Two-mode spin-squeezing as a resource for multiparameter quantum estimation 
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We study the estimation of a two-phase rotation in a finite dimensional spin system. We consider two-mode 
spin squeezed probe states and discuss the role of two-mode spin squeezing in this kind of estimation. After 
deriving the ultimate bounds posed by quantum mechanics for such probe states, we present a measurement 
scheme able to exploit two-mode spin squeezing in order to improve the measurement sensitivity. We discuss 
the case where, even in the presence of two-mode spin squeezing, the coherent state strategy yields a better 
performance and thus the standard quantum limit is not overcome. Then, by considering optimized two-mode 
squeezed probe states, we show that our measurement scheme always beats the standard quantum limit, and 
approaches the ultimate bounds for large values of the spin dimension. 
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I. INTRODUCTION 

Continuous variable (CV) squeezing EE), a ^ so known as 
bosonic squeezing, has played a major role in the develop- 
ment of modern quantum optics, and it is now safe to say 
that a good level of understanding has been reached regard- 
ing its fundamental and practical importance [3|. Its study 
has been somewhat facilitated by the fact that the uncertainty 
relation between two canonically conjugate operators is state- 
independent, hence a universal notion of squeezing arises nat- 
urally and its quantification via a 'squeezing parameter' is es- 
sentially unique (up to monotone re-scalings). 

More recently, increasing attention has been devoted to the 
finite-dimensional counterpart of CV squeezing, that is, Spin 
Squeezing HE!. This has been mainly motivated by the study 
of correlations and entanglement between particles 00, and 
by the quest to increase the precision of measurements in 
experiments involving finite-dimensional systems (3] [9]. In 
the study of spin squeezing, the first noticeable difference as 
compared to the CV case is that the uncertainty relation be- 
tween two angular momentum operators is state dependent, 
hence several inequivalent squeezing parameters can be de- 
fined (6j. Typically, the usefulness and relevance of each 
one of these quantifiers depends on the considered applica- 
tion, and the performance of different tasks will be related to 
different squeezing parameters. To mention some notable ex- 
amples, spin-squeezing and its various quantifiers have proven 
useful tools for the detection of quantum entanglement in 
many-body systems||71[8), and for increasing the precision of 
instruments such as Ramsey spectrometers ||5][9l[T0), atomic 
clocks [5, 8 1 and ultra-sensitive magnetometers ifTTIl . Even 
applications as broad as quantum computation and quantum 
simulation have been shown to benefit from the presence of 
spin squeezing [12|, and the latter is becoming progressively 
easier to be generated and measured experimentally: in par- 
ticular in atomic ensembles fl3l . Bose-Einstein condensates 
[14 1 and atomic chips fl5l . 

Considering bipartite systems, one can also introduce the con- 
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cept of two-mode spin squeezing (TMSS), as the generaliza- 
tion of its bosonic counterpart. Two mode spin squeezed states 
exhibit squeezing in their collective spin degrees of freedom, 
but not in the local ones. TMSS has been widely investigated 
lfT6l - l20 l and for pure states it has been proven to be equivalent 
(up to local unitaries) to bipartite entanglement for systems 
with equal spin lETl . 

Here, we are particularly interested in the application of spin 
squeezing to the rapidly developing field of quantum metrol- 
ogy, which studies how the peculiar features of quantum me- 
chanics affect the achievable precision in the estimation of 
one or more parameters. In doing so, we shall use the theo- 
retical tools provided by quantum estimation theory 12211231 . 
which gives fundamental bounds to the achievable accuracies 
in terms of the Quantum Fisher Information (QFI). 
In the context of quantum metrology with bosonic systems, 
squeezing is well established as being a desirable feature, able 
to enhance the performance of single-parameter estimation 
protocols 1 24 -28 1 . On the other hand two-mode squeezed 
states have been proven to be useful resources for several 
estimation problems such as the joint estimation of loss and 
temperature of a bosonic channel l29l [30l or the characteri- 
zation of both parameters defining a displacement operation 
I13T1 |32l . Similarly, spin squeezing has been identified as a 
crucial resource in several single-parameter estimation proto- 
cols, which typically boil down to determining the angle of 
rotation of a spin around some known axis Q. 
In this paper, we contribute to these developments by inves- 
tigating the resource power of TMSS in a multi-parameter 
quantum estimation problem. The estimation of multiple pa- 
rameters is receiving increasing attention in the literature (see 
for example 122] HU |30l El [3J1 ESJ [37ti42l). This kind of es- 
timation is fundamentally different from the single-parameter 
case, as different parameters may be associated to mutually 
non-commuting operators. Hence, the theoretical bounds pro- 
vided by local estimation theory may not be achievable, as 
they may assume measurement strategies that go beyond the 
operations allowed by quantum mechanics. 
More specifically, we study a broad class of two-mode spin 
squeezed states. We show how these states allow the efficient 
estimation of two unknown phases characterizing a spin ro- 
tation, by making use of a simple strategy in which the rota- 
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tion is applied to only one of the two subsystems, followed 
by a projective measurement on the local spin operators of 
both parties. By accurate numerical calculations we show that 
two-mode spin squeezing is indeed a key resource in the con- 
sidered problem, allowing to beat the standard quantum limit 
(SQL) which, for spin systems, is obtained by considering 
the estimation power of coherent spin states (CSS). Moreover, 
we show how our considered strategy becomes nearly-optimal 
when the dimensionality of the spins becomes large. The pa- 
per is organized as follows: in the next section we introduce 
the concept of spin squeezing, giving its definition in terms of 
a metrological application, along with the definition of two- 



mode spin squeezing. In Sec. Ill we provide a brief introduc- 
tion to local quantum estimation theory, both for the single- 



and the multi-parameter case. In Sec. IV we discuss in detail 
the estimation problem at the center of our investigations. We 
propose a measurement strategy for the estimation of the rel- 
evant parameters, and we compare our results to the ultimate 
limits posed by quantum mechanics. Sec. [V] concludes the 
paper with some remarks. 



II. SPIN SQUEEZING 

Let us consider a spin- j system, characterized by the angu- 
lar momentum operators {J x , J y , J z }, which satisfy the com- 
mutation rule 



(1) 



where a, /?, 7 denote the components in any three orthogonal 
directions, and e Q ^ 7 is the Levi-Civita symbol. This implies 
an uncertainty relation 



AJ X AJy > 



(2) 



where AJ 2 = (J 2 )~(J a ) 2 and (O) = Tr[gO] denotes the ex- 
pectation value of an operator on a given quantum state g. To 
better understand spin squeezing, it is useful to first introduce 
the notion of coherent spin states (CSS). These are defined as 



|C(M)) = A(M)|i,j) 

where is the eigenstate of J z with eigenvalue j and 



(3) 



R(9, <p) = exp(—i(f)Jfi) — exp[i(j)(j x sin0 — j y cos6>)] (4) 

represents a rotation of angle 4> around the axis n = 
(— sin 9, cos 9, 0) T . Different definitions of spin squeezing 
have been introduced, the most used being the ones suggested 
by Kitagawa and Ueda |4|, and by Wineland et al. (see ||6) 
for a complete review of these definitions). In the following 
we will present the details regarding the latter definition (5J, 
which is more relevant for metrological applications. Let us 
start by defining the mean-spin direction (MSD) 



(J) 



\(J)\ 



(5) 



where J = (J x , J y , J Z ) T ■ Without loss of generality, let us 
consider a state with MSD along the z-axis. This state is then 
rotated by a phase <fr around the x-axis, as described by the 
unitary operator R(—ir/2, <j>). For small values of the angle cj>, 
one can show that the phase sensitivity, obtained by measuring 
the output spin operator j° ut = R(-ir/2, ^ J y R(-ir/2, </>), 
reduces to 



A<j)- 



AJ° ut 



AJ V 



iwr>i mi 



(6) 



and is thus related to the initial variance of the operator J y . 
For a coherent spin state, the fluctuations of all the spin oper- 
ators orthogonal to the the MSD are equal to (A Jn ± ) 2 = j/2 
and then one proves that 



A0SQL = 



1 



/27' 



(7) 



which is referred to as the standard quantum limit (SQL) or 
shot-noise limit. In analogy to the CV case, a quantum state is 
spin-squeezed if it exhibits lower fluctuations in one of the di- 
rections orthogonal to the MSD; according to this definition, a 
state is said to be spin squeezed if and only if its correspond- 
ing phase sensitivity is below the shot noise limit, i.e. 



A<t> < A0 SQL . 



(8) 



One can show that this can be obtained if the fluctuations 
of the spin operator J y are below a certain threshold, i.e. 

AJ y < j/2. 

Combining the uncertainty relation Q with the inequality 
j 2 > (J 2 ) > AJ 2 , which bounds the maximum variance 
on a spin observable, one also obtains the Heisenberg limit 
(HL) 



A(j) > A0HL = 



1 

27 



(9) 



posing the ultimate limit on the estimation precision of the 
phase (f> compatible with quantum mechanics. 
Let us now move on to the definition of two-mode spin 
squeezing IT6H20ll for a bipartite spin system described by 

operators J a — (J Xa , Jy a , Jz a ) T with a = {1, 2}. Sums and 
differences of spin operators belonging to different systems 
are denoted by J a ± = J Ql ± J Q2 , and obey the uncertainty 
relations 



AJ x± AJ y± > 



>y± - 2 

A state is said to be two-mode spin squeezed if 

(AJ x ^) 2 + (AJ y+ ) 2 <\(J z+ )\. 



(10) 



(11) 



Here, the reduced fluctuations are observed in the non-local 
operators J x - and J y+ . In fact, two-mode spin squeezing is 
not only a sufficient lfl9l l20l . but also a necessary condition 
for entanglement in pure states of two subsystems with equal 
spin (except for a set of bipartite states with measure zero) 

ED. 
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III. QUANTUM ESTIMATION THEORY 



IV. TWO-PHASES QUANTUM ESTIMATION 



We will give here a brief review of quantum estimation the- 
ory l22l |23~1 l43l . Let us consider a family of quantum states 
gx, labelled by a real parameter A that we aim to estimate. The 
ultimate limit to the estimation precision of the parameter A is 
given by the quantum Cramer-Rao bound 



Var(A) > 



1 



MH{\) ' 



where M is the number of measurements performed, 

H(X) = Tr[g x Ll] 



(12) 



(13) 



is the quantum Fisher information and L\ is the Symmetric 
Logarithmic Derivative (SLD) , that is, the operator satisfying 



29 A Q\ = L\Q\ + Q\Ly 



(14) 



Note that, in principle, one can always find a quantum mea- 
surement able to attain equality in Eq. (12) , hence saturating 
the quantum Cramer-Rao bound. 

Moving on to a multiparameter scenario, let us consider a fam- 
ily of quantum states g z labelled by d different parameters 
z = {z^}, fi = 1, . . . , d. The SLD for each parameter is 
defined via 



dg z 



from which one can calculate the QFI matrix H: 



H, 



We define the 



Tr 



r(S) T (S) AS) r (S) 



(15) 



(16) 



covariance matrix elements V(z)^ lL/ — 
E[Zfj,z v ] — E[Zfj]E[z„] and consider a weight (positive defi- 
nite) matrix G. Then, the multiparameter quantum Cramer- 
Rao bounds read 



tr[GV] > ^trlGtH)- 1 ] 



(17) 



where tr[A] is the trace operation on a finite dimensional ma- 
trix A and M is the number of measurements performed. We 
observe that if we choose G = 1 we obtain the bound on the 
sum of the variances of the parameters involved, 



E Var ( 



M ~ 



1 

M' 



-trfH" 



(18) 



where we have introduced the overall multiparameter sensi- 
tivity Az. Differently from the single parameter case, the 
multiparameter bound is not always achievable, since opti- 
mal measurements for different parameters may correspond 
to non-commuting observables. We remark that other (not al- 
ways achievable) bounds to the estimation precision have been 
introduced in the multiparameter case, together with the con- 
cept of most informative bound [34 -38 1. In this manuscript 
we shall focus on the standard Cramer-rao bound described 
above. 



We consider the following estimation problem: a probe 
state g undergoes a unitary evolution 



U(<p x , 4> v ) = exp{i((j> x Jy + (j> v J x )}, 



(19) 



characterized by two different phases <j) x and cf> y which we 
aim to estimate. In the case of a physical spin, such rotation 
could be the result of the interaction between the system and 
a classical magnetic field lying in the XY plane. We shall 
work in the Heisenberg picture, characterizing the action of 
the operator in Eq. ( fT9] > via 



J° 



U{<f>x, <j> y ^JU{<t> x , (j> y ) = M{4> x ,cj> y )J (20) 



where J = (J Xl J y , J Z ) T . M((f> x , <\> y ) is the 3x3 rotation 
matrix associated to the operator U (4> Xl 4> y ): its expression is 
formally equivalent to Eq. \\9) , provided that J. x and J y are 
replaced by the corresponding 3x3 generators of the rota- 
tion group. We are interested in a regime where the values of 
the parameters <j) x and <fi y are small: in this case the matrix 
M(<p x , <fi y ) can be expanded at first order in the two phases, 
yielding 



M{cf> x ,(j) y ) a 
As a consequence one has 

J° Ut R 



1 <A 
1 -< 

-<t>x <t>y 



(21) 



1 



Jy 



^xJz 



(22) 
(23) 



showing that by measuring the angular momentum operator 
J x and J y after the rotation one gains information about the 
corresponding phases. 

Following the results obtained for multi-parameter estimation 
by using two-mode squeezed entangled states in continuous- 
variable systems |29, 30, 32), we will consider the case where 
the probe state g is a bipartite system and each subsystem 
a is described by the angular momentum operators vector 

J a = {Jx a , Jy a , Jz a ) T ■ The estimation strategy is described 
in Fig. [TJ a bipartite state is prepared and the rotation oper- 
ator U((f> x , <t> y ) is applied on one of the two arms. After the 
rotation, spin measurements are performed on the two sub- 
systems. Before specifying the measurements that have to be 
performed, let us focus on the probe state. We look for bi- 
partite states which exhibit two-mode spin-squeezing, that is, 
they satisfy Eq. ( fTTj ), and in order to do this we follow the 
approach of Ref. [ 19] . We thus consider states in the Schmidt 
form 



l*in) 



m=j 



m=-j 



$m |i, m) \j, m) 



(24) 



4 



M/2 M/2 
copies copies 




Figure 1: Measurement strategy for a two-phases measurement. A 
two-mode spin squeezed state is used as a probe. On one arm the 
two-phases rotation is performed; after the rotation measurements of 
alternatively J X1 and J X2 or J yi and J V2 are performed on the two 
output systems. The results are then combined to obtain J x - and 

Jy+. 



From Eq. ( |24) , one can immediately deduce the following re- 
lations between the expectation values of local spin observ- 
ables: 



(Jz 



(J, 



< 4) = <4 2 > = <4) = <4) = o 



(25) 
(26) 

Within the class of states defined by Eq. p4) , we can identify 
two-mode spin squeezed states by fixing a real parameter fi, 
solving the eigenvalue problem 



(27) 



and choosing the eigenvector corresponding to the minimum 
eigenvalue v min . When this procedure is repeated for differ- 
ent values of /i, it is found that the corresponding output state 
yields different values of (J z +) = 2(J Zl ). Hence, the ex- 
pectation value ( J Zl ) can be used to label the two-mode spin- 
squeezed states produced via such algorithm [19|. If at least 
two coefficients $ m are not null, it is known that the state will 
exhibit two-mode spin squeezing l2ll . Applying Eq. ( |2"0) to 
the first subsystem, together with the small angles approxima- 
tion [Eq. (21)], we obtain 



Jz 



out \ 

y+ 1 



Jz 



(28) 



As a consequence, we can obtain information about the phase 
values from the measurement of the two observables J x _ and 
J y +. Moreover, since the variances of these observables are 
not changed by such rotation, i.e. 

(AJ°^) 2 w (AJ X -) 2 , 



(AJ°f) 2 « (AJ y , 



(29) 
(30) 



one can use the TMSS property ( fTT) to perform a precise esti- 
mation of the two parameters. In particular, the intrinsic phase 
sensitivities of the probe state read 



A<j) x 



AJ X - 

1(4)1 



A(j)y = 



\{J*)\ 



(31) 
(32) 



The estimation strategy in Fig. [TJ resembles to the one used 
in [32 1 for the two-parameter displacement estimation with 
continuous-variable systems. However, in that case one is 
able to associate the two parameters to two commuting ob- 
servables (at the cost of adding some extra noise), which can 
then be measured simultaneously. Here, the two observables 
J°^ and J°+ commute only on average, and in such a case 
the question of simultaneous measurement becomes nontriv- 
ial ll33l . We then choose a measurement strategy which is 
motivated by simplicity: given M copies of the probe state, on 
M/2 of those we perform the measurement of J Xl 2 , while the 
operators J Vl 2 are measured on the remaining M/2 copies. 

Then, M/2 values for J x _ and J y + are obtained respectively 
by subtracting or summing the experimental outcomes. When 
M 3> 1, the Central Limit Theorem predicts the total experi- 
mental variance of the estimated parameters, as per 



Var(^) + Var(^) = — {A<& + A^) 



(33) 



A comparison of this expression with the Cramer-Rao bound 



in Eq. ( 18 1, which is also expressed in terms of the total num- 
ber of measurements M, leads us naturally to define the effec- 
tive phase sensitivities 



A<pf = V2A(/) X A4>f = V2A<j> y 



(34) 



Then, the performance of our estimation protocol can be quan- 
tified via the two-phase sensitivity 



A$ 



efH2 _u (A0f ) 2 



(AJ X _)2 (Ajy+f 



\ V K4)l 2 



l(^)l 2 



(35) 



Our aim is now to compare the phase sensitivity in Eq. ( 35 1 



to the ultimate bounds imposed by quantum estimation the- 
ory. The HL, representing the best sensitivity compatible with 
quantum mechanics, can be easily calculated by noting that 
the uncertainty relations in Eq. ( fTO) , together with Eq. p5| ) 
and the condition A J a ± < 2j, imply AJ^_ > \(J Zl )\/2j 
and also AJ y+ > \(J Zl )\/2j, so that the minimum value that 



(36) 



Eq. (|35|) can take compatibly with these constraints is 



A$ 



HL 



The SQL represents the estimation sensitivity which is achiev- 
able with 'classical' resources, that is, coherent spin states. 
One may be tempted to calculate such quantity by evaluat- 
ing Eq. ( |3~5j ) for a CSS such as \j,j) O \ j,j)- This, however, 
would underestimate the resource power of CSSs: we recall 
that in our protocol the spin rotation (jT9j is applied only to 
the first subsystem, so that the measurement of the second 
subsystem cannot yield information about such rotation when 
the resource state is a two mode CSS, and it just results in 
extra noise being added to the experimental data. Then, to 
obtain the correct SQL one has to simply neglect the second 
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subsystem, and consider the estimation problem of a spin-j 
CSS to which the rotation ([19} is applied, followed by the 
measurement of either J x or J y (each observable being mea- 
sured on half of the copies). Following the same steps that 
lead to the single-parameter SQL in Eq. one derives the 
two-parameter SQL as per 



Ad> 



SQL 



V2 



(37) 



0.8 
0.6 
0.4 
0.2 



<J zl > 



Finally, one can evaluate the Cramer-Rao bound for the class 



of probe states described by Eq. (24 1. The QFI matrix for the 
two-phases estimation problem reads 



H 



J 2 



(jxi Jyi ~^~Jy\ Jx\ \ /2 



7 2 



(38) 

where we have used Eq. (|26]>. As a consequence, one obtains 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



(A$) 2 > Tr[H _1 ] 



1 



P+3-{Jl) 



(39) 



<J zl > 



The Cramer-Rao bound in Eq. (39i would suggest that the 



state yielding the best two-phase sensitivity should have 
(J|j) = 0, thus implying (J 2l ) = 0. However, by apply- 
ing our measurement strategy to such a state, it is easy to see 
that no information can be gained about the two phases <fi x , 4> y 
[see Eq. <|28)], so that one obtains a divergent phase sensitiv- 
ity. This clearly implies that such limit cannot be reached. 
Hence, to find the optimal probe state for a given spin number 



j we can directly study the phase sensitivity A$ in Eq. (35 1, 



as a function of the parameter ( J Zl ) (we recall that such mean 
value uniquely labels our probe states). The typical behaviour 
of the phase sensitivity as a function of ( J 2l ) is shown in Fig. 
We observe that even if the choice ( J 2 , ) =0 would lead to 



Figure 2: (Red upper solid line) Phase sensitivity A<£> and (blue 
lower solid line) Cramer-Rao bound for the state in Eq. [24] as a 
function of (J Z1 ), along with the SQL (dashed green line) and the 
HL (dotted black line), for different values of the spin dimension j. 
Top: j — 5; bottom: j = 7. 



clearly shows that bipartite entanglement, even if it can be ex- 
ploited to the enhancement in the estimation precision, is not 
the sole resource for this particular task, an important role be- 
ing played by the detailed structure of the considered states. 

It can be observed that the considered two-mode squeezed 
states, combined with our measurement strategy, allow us to 
beat the SQL in a wide parameter range, with the only excep- 



the Cramer-Rao bound in Eq. ( 39 1. 

This behaviour is confirmed by observing Fig. [3] where 
the phase sensitivity is plotted against the entanglement of the 
probe state, quantified via the normalized von Neumann en- 
tropy of the reduced state 



E[*in] = 



S(Tr B [|* in )(* in | 



a diverging A<1>, the optimal probe state corresponds to a value 

of ( J 2l ) which is close to zero,_as it is somehow suggested by tion of the two regions ( J 2l ) ~ and (J 2l ) < j. Note in 

particular that when ( J Zl ) — » j the probe state converges to 
a CSS, while the phase sensitivity does not converge to the 
SQL: our measurement strategy is in fact sub-optimal when 
applied to CSSs (see derivation of the SQL above), and it be- 
comes advantageous only when the TMSS property compen- 
sates for the extra noise introduced by measuring the second 
subsystem. In Fig. [4[ we plot the optimized phase sensitivity 
as a function of the spin number j, comparing it to the Cramer- 
Rao bound, the Heisenberg limit and the Standard quantum 
limit. We see that our estimation protocol yields a phase sen- 
sitivity which is well below the standard quantum limit. We 
also note that our strategy appears to be closer to optimality 
as the spin number j is increased. Indeed, while for small 
values of j there is a visible gap between our phase sensitivity 
and both the Heisenberg limit and the Cramer-Rao bound, this 
difference significantly decreases upon increasing the dimen- 
sion of the Hilbert space. We recall that for multi-parameter 



log 2 (2j + l) ' 
with S(g) = -Tr[£log 2 £>]. 



(40) 
(41) 



Given the mentioned equivalence between two-mode spin 
squeezing and entanglement 11211 . and the pivotal role of the 
latter in quantum technologies, one may suspect that entangle- 
ment itself might be the essential resource in the considered 
protocol. However, the examination of Fig.[3]reveals that this 
would be an oversimplification: even though the minimum of 
A<E> is found in the region of high entanglement, the phase 
sensitivity rapidly diverges in correspondence with the max- 
imum entanglement achievable, which in fact corresponds to 
the 'phase insensitive' state ( J 2l ) = discussed above. This 



estimation the Cramer-Rao bound obtained in Eq. ( 39 1 is not 
always achievable (see Sec. Ill i. This is evident for low val- 



ues of j, where the Cramer-Rao bound is below the Heisen- 
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berg limit. Similarly, the two-parameter Heisenberg limit in 
Eq. ( |36] > is not always achievable in principle (the optimal 
states for different phases may be different and the two opti- 
mal measurements may not commute). In view of the above 
analysis, the measurement strategy we have proposed can be 
considered nearly optimal for the studied estimation problem, 
when the spin number j is large enough. 



V. CONCLUSIONS 



Figure 3: Phase sensitivity A$ (solid lines) as a function of the nor- 
malized entanglement EbJ/in], along with the SQL (dashed lines) for 
different values of the spin dimension j. From top to bottom: j — 3 
(green); j = 5 (red); j = 7 (blue). 
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Figure 4: Phase sensitivity A$ (Red circles) and corresponding 
Cramer-Rao bound (blue triangles) of optimized probe states as as a 
function of the spin dimension j, along with Heisenberg limit (black 
solid line), standard quantum limit (green dashed line). 



We have studied a multiparameter quantum estimation 
problem in which a two-phase spin rotation is characterized 
by using two-mode spin states as probes. We have presented 
a measurement strategy which, interestingly, is sub-optimal 
when the resource states are classical, while at the same time 
is able to exploit the presence of two-mode spin squeezing 
to overcome the standard quantum limit. By optimizing the 
probe states for each spin number j, we have shown how our 
estimation strategy can approach the precision limits imposed 
by quantum estimation theory, as the spin number is increased. 
Our findings can be of general interest to the field of multi- 
parameter quantum metrology with spin-squeezed resources. 
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